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Recent scanning tunneling microscopy in the superconducting regime 
of two different cuprate families has revealed unidirectional bond-centered 
modulation in the local electronic density of states. Motivated by this 
result we investigate the emergence of modulated d-wave superconductivity 
coexisting with charge domains that form along one of the crystal axes. 
While detailed stripe profiles depend on the used form of the Gutzwiller 
factors, the tendency towards a valence bond crystal remain robust. We 
also find closely related stripe phase originating from the staggered flux 
phase, a candidate for the pseudogap phase of lightly doped cuprates. 

PACS numbers: 74.72.-h, 74.20.Mn, 74.81.-g, 75.40.Mg 

1. Introduction 

It is now well established that the simplest model proposed to describe 
the physics of the high-Tc superconductors, the so-called t-J model 

W = E + h.c.) + S, • S„ (1) 

where cj^. = (1 — ni-a-)c^^^ is the Gutzwiller projected electron operator and 
riia is the particle number operator, yields apart from the true long-range 
magnetic order, characteristic of the undoped parent Mott insulators, an 
array of quantum SU(2)-invariant ground states [2]. In fact, it is quite nat- 
ural to expect that strong quantum fluctuations arising from both low spin 
S" = 1/2 of the copper ions and the two-dimensional nature of the Cu02 
planes, should lead to quantum disordered states with only short-range an- 
tiferromagnetic (AF) spin correlations. The most famous example of such 
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states is a resonating valence bond (RVB) phase [3]. Remarkably, Ander- 
son's RVB theory based on a Gutzwiller projected BCS trial wave function, 
which parameters are usually determined either by using renormalized mean 
field theory (RMFT) [3] or by Variational Monte Carlo (VMC) method [5], 
not only predicted correctly the d-wave symmetry of the superconducting 
(SC) order parameter [6], but in addition, it reproduced experimental dop- 
ing dependence of a variety of physical observables in the SC regime [7]. 

Moreover, the tendency towards valence bond amplitude maximization 
might enhance charge and spin stripe correlations in the d-wave RVB state. 
The presence of charge and long-range spin stripe order has been detected 
in neutron scaterring experiments and confirmed in resonance x-ray scatter- 
ing in a few special cuprate compounds, namely La1.6-xNdo.4Srx.CuO4 and 
La2-xBaxCu04 [8]. However, such stripe order competes with supercon- 
ductivity and thus strongly reduces Tc In contrast, recent scanning tun- 
neling microscopy (STM) on different cuprate families Ca2-xNaa;Cu02Cl2 
and Bi2Sr2Dyo.2Cao.8Cu208+5, has revealed intense spatial variations in 
asymmetry of electron tunneling currents with bias voltage that forms uni- 
directional domains coexisting with inhomogeneous d-wave superconduc- 
tivity [To]. In particular, it has been found, that the asymmetry occurs 
primarily at the oxygen sites being indicative of a short-range bond- centered 
charge pattern with a period of four lattice spacings. In this paper we 
show that the bond-centered modulation observed in the STM experiments 
might be naturally interpreted in terms of a valence bond crystal, i.e., spin- 
rotationally invariant phase with spatially varying bond charge hopping and 
a concomitant modulation of short-range AF correlations [llj . 



We begin by discussing RMFT of the t-J model applied to the case 
with homogeneous charge distribution. In this approach, the Gutzwiller 
projection removing double occupancy is handled with statistical weight 
factors gt = 2x/(l + x) and gj = 4/(1 + x)'^ which account for different 
probabilities of hopping and superexchange processes in the projected and 
unprojected wave functions. Hence the mean-field Hamiltonian reads, 



with the Bogoliubov-de Gennes self-consistency conditions for the bond- 



Xji = and pair-order Aji = {Cj_^c^^^) = (c^ _^Cj. ^) parameters in 
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Fig. 1. Electronic structure along the main directions of the Brillouin zone of: (a) 
d-wave RVB/SF (x = 0); (b) d-wave RVB and (c) SF phase (both for x = 1/8). 



the unprojected state. Hereafter, we shall assume a typical value t/J = 3. 

Even though the original proposal for the high-Tc superconductivity was 
the s-wave BCS wave function, it immediately turned out that it is the d- 
wave BCS state with Ajj = itA for the nearest-neighbor pairs along the 
X (y) axis, respectively, which gives the lowest energy [6j. At half-filling, 
such a phase is equivalent to the staggered flux (SF) state with complex 
Xij = Ixl exp[(—iy'^~^^yi(p] yielding circulating currents whose chirality al- 
ternates from plaquette to plaquette [12]. In this limit, the t-J model re- 
duces to the Heisenberg Hamiltonian with the local SU(2) gauge symmetry 
corresponding to the following particle-hole transformation. 



at 
-fit 



f3i 
a. 



(3) 



with a^a* -|- /3i/3* = 1. It mixes an |-spin particle with a j-spin hole and 
hence decouplings in terms of A or x become indeed equivalent. In order to 
appreciate this better let us consider the related Hamiltonian matrices using 
the Bogoliubov-Nambu formalism with ry^ = (c^p ^^-Vii) the d-wave RVB 



phase and r/k 



with Q = (vr, vr) for the SF phase: 



^RVB 



Ak ek + 



-Ek - /" *Xk 
-«Xk Ek - 



(4) 



where ek = (^5 + p g-^ Rex)l+, Ak = fJg-' Xk = {Jg^'lmxi- with 
7± = 2(cos /ca; lb cosky). Therefore, the corresponding spectra are given by: 



^k 



^k' 



and 



f^±\^4 + xl- (5) 



At half-filling (/x = 0), one finds A = x = 0.169 (|x| = 0.239 and (p = 7r/4) 
for the d-wave RVB (SF) phase, respectively. Hence, in the latter case 
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Fig. 2. Doping dependence of: (a) pairing amplitude A (open squares) and SC order 
parameter Age (solid squares) in the d-wave RVB phase, as well as (b) plaquette 
flux <i> (open circles) and spin correlations S (solid circles) in the SF phase. 



Rex = Imx = 0.169 and, since = Xk, both spectra become degener- 
ate. As shown in Fig. [Dj^a), the key feature of the obtained spectrum is 
that the energy gap vanishes hnearly along the S = (7r/2,7r/2) point form- 
ing a cone-like dispersion. While this cone remains pinned to the Fermi 
surface in the d-wave RVB phase, finite doping takes the node of the SF 
order away from the Fermi surface and opens hole pockets around the S 
point [see Fig. [T]||b,c)]. Nevertheless, both excitation spectra remain similar, 
which makes the SF phase an excellent candidate for the normal pseudogap 
phase that emerges below a characteristic temperature T*. Moreover, short- 
range staggered orbital current-current correlations have been found in the 
Gutzwiller-projected d-wave RVB phase [13j, in the exact ground state of 
the t-J model with a negative two-hole binding energy [T3], as well as by 
analyzing motion of a hole pair in the AF background [15j . 

Next, as shown in Fig. [5]^a), variational parameter A is the largest at 
X = and vanishes linearly with doping. In contrast, g*^ is an increasing 
function of x so that the resulting SC order parameter Asc = re- 
produces qualitatively the SC dome. Finally, Fig. EJ^b) depicts the doping 
dependence of the fictitious flux (in unit of the flux quantum) defined by a 
sum over the four bonds of the plaquette = X](ij)gn ©ij well as the 
AF spin correlations S = — fff'^lxP- Here, the appearance of a finite flux at 
a; ~ 0.15 clearly strengthens S with respect to the Fermi liquid state where 
X is entirely real. 



3. Unidirectional charge order 

We turn now to the discussion of bond-centered (with a maximum of the 
hole density spread over two-leg ladders) inhomogeneous RVB (chiral) states 
derived from the parent d-wave RVB (SF) phases, respectively. Hereafter we 
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Fig. 3. (a,b) Hole density n^i, (c,d) SC order parameter A^^, (e,f) spin correlation 
Sf , and (g,h) bond charge , found in the ttDRVB phase. Top (bottom) panels 
depict the results obtained using original (modified) Gutzwiller factors; solid (open) 
circles in panels (c-h) correspond to the x (y) direction, respectively. 



refer to the former as vr-phase domain RVB phase (vrDRVB), as it involves 
two out-of-phase SC domains (see also Ref. [IS]), separated by horizontal 
bonds with vanishing pairing amplitudes, named as "domain wall" (DW), 
where Aj-,- gains a phase shift of tt. Similarly, due to the existence of DWs 
which act as nodes for the staggered current and introduce into the SF order 
parameter a phase shift of vr, we refer to the latter as ttDSF state. 

We consider both original (q = 0) and modified (g = 1) Gutzwiller 
factors depending on local hole densities nhi, 

J ^ 4{1 - nhi){l - Hhj) 

a,, + qiSnMnhjPij (2) + 16/3+ (4)] ' ^ ^ 



^nhiUhjil - nhi){l - rihj) 




with = {l—n\^){l—n\-) and j3fj{n) = lAjjl^ztlxijl""- Note, however, that 
the Gutzwiller renormalization scheme becomes substantially more compli- 
cated in the case of inhomogeneous charge distribution as the local density 
may change before/after projection [17j. As a consequence, Eqs. ([6]) and ([7]) 
may provide only an approximate way of the Gutzwiller projection. Finally, 
using the unit cell translation symmetry p!8], calculations were carried out 
on a large 256 x 256 cluster at a low temperature fiJ = 500. 

The corresponding stripe profiles in both phases shown in Figs. [3]and[l] 
are clearly a compromise between the superexchange energy Ej and kinetic 
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Fig. 4. (a,b) Hole density Uhi, (c,d) modulated flux (e,f) spin correlation Sf , 
and (g,h) bond charge , found in the ttDSF phase. Top (bottom) panels depict 
the results obtained using original (modified) Gutzwiller factors; solid (open) circles 
in panels (e-h) correspond to the x (y) direction, respectively. 



energy Et of doped holes. On the one hand, a reduction of the SC or flux 
order parameters (the latter known to frustrate coherent hole motion ^19j ) 
enables a large bond charge hopping Tf' = 2g\ ^_^yRe{xi,i+y} along the DWs 
as in the usual stripe scenario [18]. On the other hand, it simultaneously 
results in the suppression of the AF correlations Sf = — |9/j_|_2.(|xi,i+a;P + 
along the transverse bonds. 

However, a closer inspection of Figs. [3] and [J] as well as Table I indicates 
that this competition is especially subtle in the vrDRVB state. Indeed, in- 
stead of increasing hole level in the SC areas in order to reinforce the SC 
order parameter, the system prefers a more spread out charge distribution, 
which suggests that the d-wave RVB state is less disposed to phase separa- 
tion than the SF one where in fact also other charge instabilities have been 
found (201 121j . Moreover, as listed in Table I, both Et and Ej are reduced 
with respect to the uniform d-wave RVB superconductor. In contrast, vrDSF 
phase fully optimizes both energy contributions simply by expelling holes 
from the regions between the stripes and accommodating them at the DWs. 
Indeed, a low local doping level strengthens plaquette flux which reaches 
the value <!>□ ~ 0.35 expected for x ~ 0.1 [see Fig. [2)^b)]. Moreover, the 
holes accommodated at the DWs enhance locally Gutzwiller factors g\j and 
allow the phase to retain a favorable Et. Taken together, these two effects 
are responsible for a much stronger charge modulation of the ttDSF phase 
as compared to its vrDRVB counterpart. 
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TABLE I 

RMFT kinetic energy Et, magnetic energy Ej, and free energy F of the locally 
stable phases: ttDSF, SF, ttDRVB, and d-wave RVB one at a; = 1/8. 



original g.y 


modified g.ij 


phase 


Et/J 


Ej/J 


F/J 


Et/J 


Ej/J 


F/J 


ttDSF 


-1.0252 


-0.4320 


-1.4572 


-0.8514 


-0.4269 


-1.2783 


SF 


-1.0345 


-0.4246 


-1.4591 


-0.8622 


-0.4230 


-1.2852 


ttDRVB 


-1.0160 


-0.4607 


-1.4767 


-0.8719 


-0.4518 


-1.3237 


RVB 


-1.0232 


-0.4838 


-1.5070 


-0.8863 


-0.4784 


-1.3647 



Unfortunately, the total RMFT energy in both phases differs substan- 
tiahy from the one obtained within the VMC scheme: Ed^-^/J — —1.34 
{Ebsf/J ^ -1.33) for the vrDRVB (vrDSF) phase, respectively [H]. We con- 
sider therefore the so-called modified Gutzwiller factors where the effects of 
the nearest- neighbor correlations Xij and Aij are also included [22j. First 
of all, one observes that the inclusion of the intersite correlations weakens 
(strengthens) stripe order in the ttDRVB (ttDSF) phase, respectively. In- 
deed, in both cases the holes are ejected from the regions in between stripes 
into the DWs defined as nodes of the SC/flux order parameter. Similar 
tendency towards accommodating the holes at the DWs has also been es- 
tablished in the VMC calculations [23]. Furthermore, while the short-range 
AF correlations remain either unaltered or they are changed in such a way 
that Ej remains almost constant, modified Gutzwiller factors mainly renor- 
malize bond charge hopping. As a consequence, the total energy in both 
phases approaches the one found in the VMC scheme (see Table I). 

In summary, we believe that our results provide insights into the for- 
mation of the recently found bond-centered charge order that coexists with 
modulated d-wave superconductivity. Moreover, we expect a further en- 
hancement of the proposed valence bond crystal near impurities that break 
the space group symmetry of the t-J Hamiltonian by producing a modu- 
lation in the magnitude of the superexchange coupling |23]. In fact, it has 
recently been argued that the dopant-induced spatial variation of the atomic 
levels indeed strengthens locally the AF superexchange interaction [25] . 
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